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In the context of a charge qubit under continuous monitoring by a single electron transistor, we 
propose an unraveling of the generalized quantum Markovian master equation into an ensemble 
of individual quantum trajectories for stochastic point process. A suboptimal feedback algorism 
is implemented into individual quantum trajectories to protect a desired pure state. Coherent 
oscillations of the charge qubit could be maintained in principle for an arbitrarily long time in 
case of sufficient feedback strength. The effectiveness of the feedback control is also reflected in 
the detector’s noise spectrum. The signal-to-noise ratio rises significantly with increasing feedback 
strength such that it could even exceed the Korotkov-Averin bound in quantum measurement, 
manifesting almost ideal quantum coherent oscillations of the qubit. The proposed unraveling and 
feedback protocol may open up the prospect to sustain ideal coherent oscillations of a charge qubit 
in quantum computation algorithms. 

PACS numbers: 42.50.Dv, 42.50.Lc, 05.40.-a, 73.23.-b 


I. INTRODUCTION 

The advent of quantum information technologies is cre¬ 
ating a considerable demand for strategies to manipulate 
individual quantum systems in the presence of noise 111]. 
Recently, state-of-the-art nano-fabrication has made it 
possible to accurately monitor a single quantum state in 
a continuous manner M. It opens up the prospect for 
physically implementing the measurement-based (close- 
loop) feedback control of a quantum system, in which 
the real-time information of the detector’s output is ex¬ 
tracted and instantaneously fed appropriate corrections 
back onto it HI- So far, a variety quantum feed¬ 
back protocols have been proposed [lfli (HI • These con¬ 
trol strategies turn out to be efficient for potential ap¬ 
plications ranging from the realization of a mesoscopic 
Maxwell’s demon fTjl - ITjf and purification of a charge 
qubit iTHt - TTr^ , to the freezing of charge distribution [l8[ 
in full counting statistics (FCS) [H, [2(| . 

In contrast to the classical feedback, by which one 
might acquire as much system information as possible, 
in quantum feedback control the measured system is in¬ 
evitably disturbed in an unpredictable manner, and the 
amount of information is fundamentally limited by the 
Heisenberg’s uncertainty principle. The essence of the 
involving measurement process is the trade-off between 
the acquisition of information and the backaction induced 
dephasing of the quantum state f2Tl - [2d] | . It is therefore 
suggested to consider a continuous weak quantum mea¬ 
surement, in which information gain and disturbance of 
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the measured system could reach a balance. The Lind- 
blad master equation is, in general, utilized for contin¬ 
uous measurement and quantum feedback control. The 
unique advantage is its clear physical implications and 
simplicity to be unraveled into individual quantum tra¬ 
jectories, in which quantum feedback operations might 
be appropriately implemented M. 

It was revealed, however, the Lindblad maste r eq uation 
may not satisfy the detailed balance relation |2f|. Fur¬ 
thermore, the dynamics of a quantum system does not 
necessarily obey the Lindblad master equation in some 
parameter regimes (for instance, the applied voltage is 
comparable to the internal energy scales of the quan¬ 
tum system). A generalized quantum Markovian mas¬ 
ter equation (G-QMME) has been employed to study 
the unique features involved in various mesoscopic struc¬ 
tures, such as super-Poissonian noise in transport dou¬ 
ble quantum dot systems |27l - (29| . additional dephasing 
in an Aharonov-Bohm interferometer [30, f3l[ . as well 
as “an exchange magnetic field” in a quantum dot spin 
valve structure [32el34j . Yet, direct application of the G- 
QMME to quantum feedback control is restricted due to 
its difficulty to be unraveled into individual trajectories. 
An essential issue, therefore, is to find a recipe which is 
capable of unraveling the G-QMME with specific physi¬ 
cal meanings such that it can be applied in the realm of 
quantum feedback control. 

In this work, we propose a physical unraveling of the 
G-QMME in the context of qubit measurement by a sin¬ 
gle electron transistor (SET). It is then implemented to 
stabilize coherent evolution of the measured qubit based 
on a suboptimal feedback algorithm [35, HU. Actually, 
an unraveling scheme of a G-QMME was recently con¬ 
structed for a diffusion process [371 ] . which, however, lim- 
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FIG. 1: Schematic setup of a solid-state charge 
qubit under continuous monitoring by an SET 
detector. The charge qubit is represented as an 
extra electron confined in a double quantum dot. 
The SET detector is a single QD sandwiched be¬ 
tween the left and right electrodes, where the 
quantum dot is capacitively coupled to the qubit 
via Coulomb repulsion. Possible electron config¬ 
urations of the reduced quantum system (qubit 
plus SET QD) are displayed in (a)-(d), respec¬ 
tively. 


ited its applications only to diffusive detectors, like a 
quantum point contact (QPC). For an SET detector, the 
measurement current exhibits a typical stochastic point 
process, which comprises unique advantages of high sen¬ 
sitivity and triggering feedback operations straightfor¬ 
wardly. An unraveling scheme tailored to this process 
is therefore of crucially importance for physically imple¬ 
ment the measurement-based feedback control in the con¬ 
text of an SET detector. 

The remainder of the paper is organized as follows. We 
start with a description of the measurement setup in sec¬ 
tion El The G-QMME for the dynamics of the reduced 
system is presented in section HITl It is then followed in 
section IIVI by the unraveling of the G-QMME into an 
ensemble of individual quantum trajectories, in which 
the dynamics of the reduced system in every single re¬ 
alization of continuous measurements is conditioned on 
stochastic point process. Quantum trajectories and cor¬ 
responding stochastic tunneling events are numerically 
simulated in section [V] based on the Monte Carlo scheme. 
Section |Vl] is devoted to the implementation of feedback 
control protocol to protect coherent oscillation of the 
qubit states, which leads to the violation of the Korotkov- 
Averin bound for a sufficient feedback strength. Finally, 
we conclude in section ivm 


II. MODEL DESCRIPTION 

The setup of charge qubit under continuous measure¬ 
ment by an SET detector is schematically shown in Fig. [I] 
The charge qubit is represented as an extra electron in 
a double quantum dot. Whenever the electron occupies 
the lower (upper) quantum dot [cf. Fig.[l], the qubit is 
said to be in the logical state |a) (|/3)). The measure¬ 
ment SET is a device with a single quantum dot (QD) 
sandwiched between the left and right electrodes, where 
the QD is capacitively coupled to the qubit via Coulomb 
interaction. The energy level of the single QD is suscep¬ 


tible to changes in the nearby electrostatic environment, 
and thus can be used to sense the location of the extra 
electron in qubit. The Hamiltonian of the entire system 
reads 

H = H s + H b + H'. (1) 

The first part 

Hg = —e <j z + Q a x + (Eg + 5U)d^d, (2) 

describes the charge qubit, SET QD and their coupling, 
where d) ( d ) is the creation (annihilation) operator for 
an electron in the QD, e is the qubit level mismatch, Q 
is the qubit interdot tunneling, and a z = |ct)(oj — \/3}(/3\ 
and a x = |a}(/?| + | or) (/3| are pseudo-spin operators for 
the qubit. We assume there is only one energy level Eg 
within the bias regime defined by the Fermi energies in 
the left and right electrodes of the SET. Whenever the 
qubit occupies the logical state |/3), the SET QD “feels” 
the existence of a strong Coulomb repulsion SU (as in¬ 
ferred in Fig.[T]), which influences transport through the 
SET. It is right this mechanism that the qubit informa¬ 
tion may be acquired via the readout of the SET. 

The left and right electrodes of the SET are modeled as 
reservoirs of noninteracting electrons, with Hamiltonian 

h b= ^2 ^2 £ikc lk c ^- ( 3 ) 

fcL.R k 

Here c\ k (c^k) stands for creation (annihilation) of an 
electron in the left (£=L) or right (£=R) electrode with 
momentum k. The reservoirs are characterized by the 
Fermi functions /l/r(w), in which the involving Fermi 
energies determine the applied bias voltage, i.e. V = 
Ml — MR- Throughout this work, we set h = e = 1 for 
the Planck constant and electron charge, unless stated 
otherwise. 
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The last component of Eq. m depicts the electron tun¬ 
neling through the SET 

H ' = Y,^ kc lk d + h - c -) = E(^+ h -c-), ( 4 ) 

i,k e 

where Ft and F'j are implicitly defined. The tunnel- 
coupling strength to the electrodes is given by the in¬ 
trinsic tunneling width T^w) = 27r^ fc \ttk\ 2 5(co — £tk)- 
In what follows we assume wide-band limit, which leads 
to energy independent tunnel-couplings T L / R . 

III. GENERALIZED QUANTUM MARKOVIAN 
MASTER EQUATION 


causes prominent spin precession behavior [32l - l34j . Even 
in a double quantum dot transport system, the dispersion 
leads to bunching of tunneling events, and eventually to 
the strong super-Poissonian shot noise [H, [29[ . 

Owing to the presence of the dispersion functions, 
the G-QMME as shown in Eq. (0 is apparent a non- 
Lindblad-type master equation. How to unravel the G- 
QMME in an appropriate way such that the quantum 
feedback may be implemented in every individual quan¬ 
tum trajectories become an essential issue in quantum 
control. In the context of continuous qubit measurement 
by an SET detector, we will present a physical unrav¬ 
eling of the G-QMME [Eq. 0] exclusively for the point 
process in the following section. 


Under the second-order Born-Markov approximation, 
the dynamics of the reduced system (qubit plus SET QD) 
can be described by the following G-QMME pisj - Pill 

p{t) = -LC 0 Pit) ~ Kp(t), (5a) 

Kp{t) = ^{{d\ A { ~ ] p(t) + p(t)A {+) ]} + h.c., (5b) 

where £ 0 (‘ • •) = [Hs, (•••)] denotes the Liouvillian asso¬ 
ciated with the system Hamiltonian Eq. 0, and A.W = 

with 4 ±5 = [C^iFC) +iD ( e ± \±C)]d. 

The involving bath spectral functions C\^ (±£) are re¬ 
lated to electron tunneling through the left and right 
junctions of the SET detector 

/ OO 

d tcf\t)e ±iLt . (6) 

-oo 


IV. UNRAVELING OF THE G-QMME 

The G-QMME [Eq. 0] only depicts the dynamics of 
the reduced system under continuous measurement by 
an SET detector. In order to achieve a description of the 
readout from the SET detector, we shall resolve the num¬ 
ber of electrons transport through the SET in Eq. 0, 
and arrive at a number-resolved quantum master equa¬ 
tion [ril. Dot 

P^’Nk) = _ {i£ 0 + Uo + nt + /9 (jVLlJVa) , (9) 

where 

TZ oP ( n M = l{d^A^ P ^ N ^+p^ N ^A^d^}+h.c. 

(10a) 

represents continuous evolution of the reduced system, 
while 


Here the bath correlation functions are given by 

C ( t\t) = {F}{i)F l ) B , (7a) 

= (F e (t)F}) B , (7b) 

with (• • -)b = tre[(- • • ) P b], and P b the local thermal equi¬ 
librium state of the SET electrodes. With the knowledge 
of the bath spectral functions, the dispersion functions 
D^\±C) can be evaluated via the Kramers-Kronig re¬ 
lation 


D ( f\±C) 



7 r 



dcj 


c^\±u) 

C — LO 


( 8 ) 


where V stands for the Cauchy principal value. The dis¬ 
persion functions are responsible for the renormalization 
of the internal energy scales, which might have impor¬ 
tant roles to play in mesoscopic transport. Actually, a 
number of unique feature have already been revealed due 
to the dispersion functions. For instance, a bias depen¬ 
dent phase shift and additional dephasing of the quantum 
states in a mesoscopic Aharonov-Bohm interferometer is 
revealed [30, [3lj| . In a QD spin valve, the dispersion 
gives rise to an effective exchange magnetic field, which 


nlp {N ^ Nn) = -^{d f p (JVL+1,JVll) A[ +) } + h.c. (10b) 

7 Z£ P ^ Ntj ’ Nr ' ) = — + h.c. (10c) 

and 

n+ P {N]l ’ Nr) = -i{d t p (iVL ’ JVR+ 1 ) ^k +) } + h -C. ( 10d ) 

^rP (JV l,JVr) = + h.c. (lOe) 

describes electron transport through the left and right 
junctions of the SET detector. Here N( denotes the num¬ 
ber of electrons tunneled through the left (£=L) or right 
(£=R) junction of the SET up to time t. This equation 
cannot be solved directly, since it actually corresponds 
to an infinite number of coupled equations. We thus 
perform a discrete Fourier transformation p(xl, Xr, t) = 
E^N R el{N ^ +NRXR) P iNM (t) and finally arrives at 
a %-resolved master equation 

§- t P = jhCo + K 0 + + e ±{ ^ni) j p 


= AXUXR )p, 


( 11 ) 
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where C{\ l,Xr) is defined implicitly, and Xl/r is the 
so-called counting field in FCS [H!, Hcj. The solution to 
Eq. (1111) reads formally 

P(x L,XR,t) = e c ^ Xh ’ XR)it ~ to) p(xr,X R,to)- (12) 

We assume that electron counting begins at to such that 
p( N ^> N ^){t 0 ) = p{t 0 )6 Nhfi S NR ,o and thus p(xl,Xr,*o) = 
p(to)- By performing the inverse Fourier transformation, 
one readily obtains 

p(N L ,N R )^ _ f dXLdXR e £( XLlXR )(t-f 0 )-iy;, Nextpfo'j 

J — 7T (2 7r )“ 

= U(N h ,N R ,t-t 0 )p(t 0 ), (13) 

where the implicitly defined G(Al, A R , 6t) actually is the 
number-resolved propagator for the reduced system. An 
important figure of merit of this propagator is that it 
solely depends on the dynamic structure of the quantum 
master equation, rather than the initial reduced quantum 
state, which makes it very efficient in unraveling of the 
G-QMME into individual quantum trajectories. 

Specifically, let us consider the evolution of the reduced 
state during an infinitesimal time interval [tk,tk + d2]. 
In general, Eq. (1131) is valid for an arbitrary number of 
electrons (Al or A R ) transmitted trough the SET, de¬ 
pending on the time interval (df). Yet, here we are inter¬ 
ested in the regime where single electron tunneling events 
dominates. We take d t -C min^^ 1 ,T^ 1 ), such that the 
probability of having (Al, A r ) > 2 is negligible. Further¬ 
more, in stead of using Al and A R directly, we introduce 
two stochastic point variables dNj,(t) and dA R (f) (with 
values either 0 or 1) to represent, respectively, numbers 
of electron tunneled through the left and right junctions 
of the SET during the time interval dt. According to 
Eq. (fl3l) . given the condition of a state p(t k ) at t k , the 
state at tk + dt reads 

p (diV L) diVa)( 4 + dt) = W (dA L ,dA R ,df)p(f fe ), (14) 

Here tunneling through the left and right junctions could 
not take place at the same time, due to the sequential 
tunneling picture of the G-QMME (JSJ). The stochastic 
point variables dNj,(t) and dA R (t) thus cannot take the 
value 1 simultaneously. 

If the measurement records are completely ignored (av¬ 
eraged over), the ensemble averaged quantum state is 
given by 

P (t k +dt)= y P idNLANR \t k +dt) 

dJV L ,dAf R 

= Y Pr(dA L ,dA R )p c (t fc + dt), (15) 

dJV L ,dAf R 

where Pr(dA R , dA R ) = tr{p( dJVL ’ dJVR )(t fe +dt)} represents 
the probability of having dA R electron tunneled through 
the left junction and dA R electron through the right one 
at the time t k + dt, and tr{- • •} denotes the trace over 
the degrees of freedom of the reduced quantum states 


(qubit plus SET QD). p c (t k + dt) = p^’ dNR \t k + 
dt)/Pr(dAL, dA R ) is the normalized state conditioned on 
the definite measurement result at t k + dt. 

In fact, Eq. (fl6l) means that if one generates dA R and 
dA R stochastically for each time interval [t k , t k + dt] and 
then collapses onto a specific state p c (t k + dt) at the 
end of each time interval, one has actually achieved a 
particular single realization of continuous measurements 
conditioned on the specific measurement results. The 
stochastic point variables dA R and dA R for single elec¬ 
tron tunneling events, respectively, satisfy 

E[dAL,(t)] = rLtr{d^p c (i)d}dt, (16a) 

E[dA R (t)] = r R tr{ctp c (t)d t }dt, (16b) 

where E[---] stands for an ensemble average of a large 
number of quantum trajectories. 

Now, it is apparent that electron tunneling conditions 
future evolution of the reduced state [Eq. (|T4j) ]. while 
real-time quantum state conditions the observed tunnel¬ 
ing events through the left and right junctions [Eq. (Il6|) ]. 
We have now successfully unraveled the G-QMME into 
an ensemble of individual quantum trajectories. Within 
this unraveling scheme, one is able to propagate the con¬ 
ditioned quantum state [/0 C (Q] and the observed result 
[dA L/R (f)] in a self-consistent manner. Furthermore, the 
unraveling of the G-QMME enables us to calculate the 
noise power spectrum via the stochastic formalism. The 
fluctuations in the detecting current are characterized by 
the two-time correlation function 

G(t) = {E[^(t+r)^(i)]-E[^(f-|-r)]E[^(Q]}| t _ ) . 00 , (17) 

where ig(t) = dNp/dt is the current through the junction 
£ ={L, R} for point process. The noise power spectrum 
of the current is then simply given by 

/ OO 

dre itJT G(r). (18) 

-OO 


V. MONTE CARLO SIMULATION FOR 
CONDITIONAL EVOLUTION 

The possible four electron configurations of the re¬ 
duced quantum system (qubit plus SET) is shown 
schematically in Fig.[l] The eigenenergies for a sym¬ 
metric qybit (e = 0) can be readily obtained as A+ ± O 
and A_ ± 0, where A± = \{8U ± \J8U 2 + Q 2 ). Here, 
we consider the case of 5U 12. Then the eigenener¬ 
gies may be greatly simplified, i.e. A + ± 12 ps A_|_ and 
A_±!l« A_. Under the application of an appropriate 
bias voltage A+ > /il > A_ > /i R , the involving Fermi 
functions in the tunneling rates are approximated by ei¬ 
ther one or zero, and the resultant quantum master equa¬ 
tion describing the reduced dynamics can be remarkably 
simplified. Let us denote the density matrix elements by 
Pjj> where j,j' = {a, b,c,d} corresponds to one of the 
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FIG. 2: Typical quantum trajectories and corresponding detection records for T2/r = 1.0 (a)-(f) and ST/T = 0.1 (g)-(l), where 
T = Tl + Tr is supposed to be a constant, and is used as the unit of energy. The initial condition is an empty SET QD with the 
qubit in the logical state a as shown in Fig. Oja), he., p c (t = 0) = |a)(a|. The measurement voltage is V/T = 5 (/Ul/r = ±2.5r), 
such that A+ > /Ul > A_ > fiR. The involving Fermi functions can be well approximated by either one or zero. The time step 
used is At=0.0ir -1 . Other plotting parameters are: e = 0, P-qT /T = 1.0, Tl = Tr = T/2, and SU/T = 10. 


electron configurations as shown in Fig.[TJa)-(d), respec¬ 
tively. The x-resolved quantum master equation m 
reads 

d _ 

~7^Paa = i^(Pab Pba) ^LPaa T 6 ^TrPcc, (19a) 

=iO(p ba - p ab ) + (e iXL r L + e iXR r R )p dd , (19b) 
d _ 

~7^Pcc = l^(Pcd Pdc ) "b & ^ L TLPaa TrPcc? (19c) 
d 

T^Pdd = i^(pdc Pcd ) (^L "b TR)Pdd 5 (I9d) 

d „ „ . „ i 

T^P&b = iO(p a a Pbb) ”b iA(p a b Pcd) ’^'hPah 

+ Qe iXL r L + p cd , (19e) 

_ _ _ i _. 

~Qj_Pcd = i^(Pcc Pdd) lA(p ab Pcd ) T ^l-A/lab 


+ ( — — r R + i SU ) Pcd ) 


(191) 


where A, the renormalization of qubit level mismatch 
arising from bath spectral functions [cf. Eq. (|8j)]. is given 
by m 


E ES* 


^=L,R ± 


27r 


1 A±-« 

2 2tt k B T 


( 20 ) 


Here U/ is the digamma function, p,t is the chemical po¬ 
tential of left (£ =L) or right (£=R) SET electrode, ks 
is the Boltzmann constant, and T is the temperature. It 
has been revealed that this renormalization has an im¬ 
portant impact on the signal-to-noise ratio in continuous 
measurement of a chage qubit. Apparently, this effect is 
absent in the Lindblad master equation. 

We are now in a position to employ the theory devel¬ 
oped in section IIVI to unravel the quantum master equa¬ 
tion (1191) into individual quantum trajectories, which will 
be used to implement feedback control of the charge qubit 
later. The numerical results for conditional state evolu¬ 
tion and corresponding tunneling events are displayed in 
Fig. [5] for different values of qubit interdot coupling f2. 
The initial condition corresponds to the charge config¬ 
uration as shown in Fig. Oja), i.e. p c (t = 0) = |a)(a|. 
For fi/T = 1.0, the qubit exhibits oscillations between 
the quantum states |a) and |b), which is stochastically 
interrupted whenever one electron tunnels into the SET 
[see Fig. m e) for dA r b[- The electron may dwell on the 
SET for a random amount of time, where it experiences 
fast oscillations with frequency ~ A+, before it tunnels 
out of the SET QD. The conditional state evolution and 
corresponding tunneling events become quite different in 
the case of suppressed interdot coupling fl/T = 0.1. One 
observes very slow oscillations between |a) and |b), as 
shown in Fig.[2Kg)-(h). Unambiguously, we find bunch- 
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FIG. 3: Unconditional measurement dynamics obtained by using the G-QMME (0 (solid curves), ensemble average over 20 
thousand quantum trajectories similar to those in Fig. [2] (symbols). For comparison, the results using Lindblad master equation 
are also plotted by the dashed curves. The initial condition is p(t = 0) = |a)(a|. The total tunneling width T = Fl + Tr is 
kept constant, and used as the unit of energy. Other plotting parameters used are: e = 0, 12/r = 0.5, V/T = 5, k^T /F = 1, 
Fl = Tr = r/2, and 5U/T - 10. 


ing of electron tunneling events through the SET; see 
Fig.0k)-(1). When the qubit stays in the logical state 
|/3), it will prevent electrons tunnel through the SET due 
to the strong Coulomb repulsion between the qubit and 
SET QD (SU/T 1). Electrons can only flow in short 
time windows where the qubit relaxes to the state |a), 
leading eventually to the bunching of tunneling events. 

The occurrence of the bunching of tunneling events is 
also manifested in the FCS of SET current. By using the 
number-resolved quantum master equation ([9]), the first 
cumulant (current) is given by 


(I) 


r L r R 
r L + 2r R ’ 


( 21 ) 


which is apparently independent of the dynamical struc¬ 
ture of the qubit. The second cumulant (I 2 ) is associ¬ 
ated to the shot noise. The corresponding Fano factor 
F = (I 2 )/(I) is given by 


F = 


r L + 4I R 

(T l + 2Tr) 2 


r L r R 


2(r L + 2r R ) 2 n 2 


O 



, ( 22 ) 


which unambiguously depends on the the qubit param¬ 
eter (f1). Here, the first term is simply the shot noise 
of a double-barrier system jjjl |46j, HH and is definitely 
below the Poisson value. The influence of qubit dynam¬ 
ics on the SET transport is characterized by the second 
term. The noise may be remarkably enhanced as the 
qubit interdot coupling H decreases, leading to bunching 
of tunneling events as shown in Fig.[2jk)- (1). Our result 
thus shows unambiguously that the intrinsic dynamics of 
the qubit may serve as an essential mechanism that leads 
to bunching of tunneling events, and finally to a strong 
super-Poissonian noise in transport through a single QD 
SET. 


To further verify the validity of our Monte Carlo 
method for conditional evolution of the reduced state, it 
is instructive to make ensemble average over a large num¬ 
ber of quantum trajectories analogous to those in Fig.O 
and then compare with the results obtained using the 
G-QMME ©. Yet, the G-QMME or the quantum tra¬ 
jectory describes the propagation of the combined (qubit 
plus SET QD) system. One thus has to trace over the 
degrees of freedom of the QD to obtain the dynamics of 
the qubit alone 


g(t) = trsET-MQ}, (23) 

where trsET{- ■ • } stands for the trace over the degrees 
of freedom of the SET QD. One may expect perform 
the trace directly upon Eq. (0; Yet, a closed form of the 
master equation for the qubit alone cannot be obtained 
without further approximations. One possible method is 
to assume extremely asymmetric tunnel junctions for the 
SET and then adiabatically eliminate the degrees of free¬ 
dom of the SET island to arrive at the reduced density 
matrix for the qubit alone (49|. Actually, such an as¬ 
sumption is equivalent to treating the SET effectively as 
a single junction detector, analogous to a quantum point 
contact. 

Here, we first numerically propagate the combined 
(qubit plus SET) system, and then trace over the de¬ 
grees of freedom of the SET QD, i.e. Eq. (l23l) to get 
the dynamics of the qubit alone. The numerical result 
is shown in Fig. [3] The reduced dynamics by using G- 
QMME (solid curves) and by ensemble average over 20 
thousand quantum trajectories (symbols) show striking 
agreement, and thus verifies the validity of the unraveling 
scheme proposed in section m For comparison, the re¬ 
sult obtained using the Lindblad master equation is also 
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FIG. 4: The dynamics of the charge qubit and corresponding degree of synchronization for different values of feedback strength 
v = 0.01 (solid curves), v = 0.1 (dashed curves), and v = 1.0 (dotted curves). The results are obtained by ensemble average 
over 20 thousand quantum trajectories. The total tunneling width T = Fl + Fr, is constant and used as the unit of energy. 
Other plotting parameters used are the same as those in Fig. [3] 


plotted by the dashed curves in Fig. [3] Notable difference 
is observed, which arises actually from the renormaliza¬ 
tion of the level mismatch A in Eq. (EU1) . 


VI. IMPLEMENTATION OF FEEDBACK 
CONTROL 

So far, we have presented a Monte Carlo method 
to simulate the continuous quantum measurement of a 
charge qubit by an SET detector. Now we are in a posi¬ 
tion to implement the feedback control of coherent evolu¬ 
tion of the charge qubit by utilizing the unraveling of the 
underlying measurement dynamics governed by the G- 
QMME, which is not necessarily of non-Lindblad form. 
In the dot state of the qubit (|a) and 1/3}), the desired 
(target) pure state under protect is 

IV’d (t)) = cos(fit)ja) + isin(f2t)|/3) (24) 

where f l is the inherent qubit interdot coupling. 

The basic idea of the feedback is to convert the detec¬ 
tor’s output into the evolution of a qubit state g c (t ), with 
g c (t) = trsET{p c (t)}. The real-time state g c (t) is then 
compared with the desired state g d (t) = |^d(i))(V , d(i)| > 
their difference is utilized to design the feedback Hamil¬ 
tonian such that their difference can be reduced in the 
next evolution step. In each successive step, the feed¬ 
back Hamiltonian acts only for an infinitesimal time in¬ 
terval d t. According to the suboptimal algorithm T>j . 
state propagation in each infinitesimal time step maxi¬ 
mizes the fidelity of g c (t) with g d (t). Specifically, let us 
consider state evolution concerning the feedback Hamil¬ 


tonian, the state g c (t + df) is given by 
g c (t + d t) =g c (t) - i[i7fb, £> c (f)]dt 

-±[H fb ,[H a> , e c mm 2 +" -, ( 25 ) 

where the feedback Hamiltonian Hr, is to be determined 
via corresponding restrictions. 

The fidelity of the state g c (t + dt) with the target state 
then is simply given by 

Fid ={ipd\ g c (t + dt)\i/j d ) 

=(iJjd\g c (t)\i/j d ) - i(^>d|[-fffd, £ ,c (f)]|f/>d)dt 

- 2<^d|[-fffb, [-Hfb, £> c (t)]]|V>d)(dt) 2 H-. (26) 

To optimize the fidelity, one should maximize the dom¬ 
inant term, i.e. the one proportional to dt. This im¬ 
poses certain constraints, e.g. the sum of the squares 
of the eigenvalues, the sum of the norms of the eigenval¬ 
ues, or the restriction on the maximum eigenvalue of F/fb, 
etc. Actually, these constraints originate from the lim¬ 
itation on the feedback strength or finite Hamiltonian 
resources. By employing the first type of constraint, 
i.e. tr q b{FF(j(} < v with tr q b{- • •} the trace over the 
qubit states, the feedback Hamiltonian is constructed as 

sm 

Hfb = g c {t)], (27) 

where n = with p = (^ d |[p c (t)] 2 |V’d} and q = 

m\g c md)?. 

To make it much more easily accessible for experimen¬ 
talists to implement the above feedback control protocol, 
we insert Eq. (l24ll into Eq. (1271) . and eventually arrive at 
a translation of the feedback Hamiltonian as 

Hfb = Kg, 


(28) 
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FIG. 5: Correlation functions of transport current through the left junction and the corresponding noise spectrum for different 
values of feedback strength v = 0.01 (solid curves), v = 0.1 (dashed curves), and v = 1.0 (dotted curves). The correlations are 
obtained by ensemble average over 8,000 quantum trajectories. The total tunneling width r = rL + FR,is set as the unit of 
energy. Other plotting parameters are the same as those in Fig. [3] 


with 


K = 


+a/?, A0<O 
-y/f, A0>O- 


(29) 


It resembles a simple bang-bang control, where the feed¬ 
back parameter between two values is simply determined 
by the phase “error” defined as A0 = <fi(t) — (f> o, with 
< t>(t) = arctan(2Im{ Q c a p{t)}/[g c aa (t) - Q c pp {t)]) the rela¬ 
tive phase between the logical states “|a)” and “|/3)”, 
and <f >o = 2f2t(mod 27 t). It is worthwhile to mention that 
the two versions of the feedback Hamiltonian protocols, 
i.e. Eq. (1271) and Eq. (1751) are totally equivalent to each. 
Yet, second one actually provided a much easier way for 
experimentalists to implement. 

In Fig.® the effect of feedback control is plotted for 
various values of feedback strength v. The feedback 
is implemented into every single quantum trajectories, 
analogous to those in Fig.® Eventually the propaga¬ 
tion of the state is obtained by an ensemble average over 
20 thousand trajectories. For sufficient large feedback 
strength (see the dotted curves for v= 1), coherent oscil¬ 
lation of the charge quit can be maintained, in principle, 
for an arbitrarily long time. To quantitatively charac¬ 
terize how close to the desired state the protected state 
reach, we introduce the synchronization degree, defined 
as D = 2tr q b{p c gd} — 1- For D = 1, it means com¬ 
plete synchronization, indicating that the state is per¬ 
fectly protected via the feedback control. The results for 
various feedback strengths are shown in Fig.HJd). Ap¬ 
parently, for a large feedback strength (v = 1), the syn¬ 
chronization degree can reach almost the maximum value 
of 1, showing thus the high effectiveness of our feedback 
scheme. 

Another important figure of merit of the coherent os¬ 
cillations is the signal-to-noise ratio of the detector’s out¬ 


put. To investigate this important feature of feedback ef¬ 
fect, we first evaluate the correlation function of current 
transport through the left junction of the SET (the right 
one gives similar results) according to Eq. (117[) . The nu¬ 
merical results obtained by ensemble average over 8,000 
quantum trajectories are displayed in Fig.® for differ¬ 
ent feedback strength v. As v increases, the correlation 
function demonstrates clearly coherent oscillation behav¬ 
ior; see Fig.®a) for v = 1.0. With the knowledge of 
the current correlation function, the corresponding noise 
spectrum can be obtained directly by performing the 
Fourier transformation, as shown in Eq. (1181) . The nu¬ 
merical results are shown in Fig.Ojd) . An enhancement 
of the feedback strength leads to a rising and sharping 
of the coherent peak located at frequency u> = 20. Ac¬ 
tually, the height of this peak relative to its pedestal is 
a measure of the signa-to-noise ratio. It was argued in 
that the signal-to-noise ratio is limited at 4, known as 
the Korotkov-Averin bound ©EH- It was revealed the 
signal-to-noise ratio of an SET detector can not reach 
the limit of an ideal detector [j52|. However, we find the 
signal-to-noise ratio can be remarkably increased in the 
case of strong feedback strength (see the dotted curve), 
even violating Korotkov-Averin bound. Our feedback 
protocol thus serves an effective method to improve the 
signal-to-noise ratio that can exceed the up bound limit 
of 4, indicating ideal qubit coherent oscillations under 
continuous measurement by an SET detector. 

Recently, a closed-loop feedback scheme was proposed 
for the stabilization of a pure qubit state by employing, 
analogously, a capacitively coupled SET. There, the real¬ 
time tunneling events through SET detector was directly 
back-coupled into qubit parameters. This purification 
process was found to be independent of the initial state 
and could be accomplished simply after a few electron 
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jumps through the SET detector. The qubit states is 
stabilized above a critical detector-qubit coupling. In 
comparison, the advantage of our feedback scheme is that 
a coherent oscillations of the qubit could be maintained 
for arbitrary qubit parameters, as long as the feedback 
strength is sufficient strong. 

VII. SUMMARY 

In the context of a charge qubit under continuous 
monitoring by a single electron transistor, we have pro¬ 
posed an unraveling scheme of the G-QMME into an en¬ 
semble of individual quantum trajectories specifically for 
stochastic point process. A suboptimal feedback algorism 
is implemented into every single quantum trajectory to 
protect a desired pure state. It is demonstrated that co¬ 
herent oscillation of the charge qubit can be maintained 
for an arbitrarily long time in the case of sufficient feed¬ 
back strength. The corresponding synchronization de¬ 


gree can reach almost close to the maximum value of 1. 
It is also observed that the signal-to-noise ratio increases 
with rising feedback strength, and remarkably could even 
exceed the well-known Korotkov-Averin bound in quan¬ 
tum measurement, reflecting actually ideal quantum co¬ 
herent oscillations of the qubit. The proposed unraveling 
and feedback scheme thus may serve as an essential mech¬ 
anism to sustain ideal quantum coherent oscillations in a 
very transparent and straightforward manner. It is also 
highly expect that it may have important applications in 
the field of solid-state in quantum computational algo¬ 
rithms. 
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